Non-exponential relaxation for anomalous diffusion 
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We study the relaxation process in normal and anomalous diffusion regimes for systems described 
by a generalized Langevin equation (GLE). We demonstrate the existence of a very general correla- 
tion function which describes the relaxation phenomena. Such function is even; therefore, it cannot 
be an exponential or a stretched exponential. However, for a proper choice of the parameters, those 
functions can be reproduced within certain intervals with good precision. We also show the passage 
from the non-Markovian to the Markovian behaviour in the normal diffusion regime. For times 
longer than the relaxation time, the correlation function for anomalous diffusion becomes a power 
law for broadband noise. 



G ; I. INTRODUCTION 

The research on the striking universality properties of slow relaxation dynamics has been driving great efforts 
in the last decades. A large and growing literature can be found in which non-exponential behaviour has been 
observed for correlation functions: supercooled colloidal systems 0, glasses and granular material 0,0], liquid crystal 
polymers 0, , hydrated proteins [|J , growth , plasmas Q and disordered vortex lattices in superconductors 
are just a few examples. Those systems present features which are similar to those found in systems with anomalous 
diffusion. The attempt to obtain response functions which are able to explain such relaxation processes is a subject 
more than a hundred years old. Rudolph Kohlrausch used stretched exponentials R(t) ~ exp[— (f/r)"] with < < 1 
to describe charge relaxation in a Leyden gas [ifj|. Later on, his son, Friedrich Kohlrausch [llj observed two distinct 
universalities: the stretched exponential with < (3 < 2, and the power law behaviour. The former behaviour is now 
known as the Kohlrausch- Williams- Watts (KWW) stretched exponential. 

There are now two main methods to describe those relaxations: a fractional derivatives approach |T3 . ll3l , ll4j . and a 
recurrence method applied to the Mori equation |l5j. However, new methods appear every day and we should mention 
£SJ ■ here the attempt of Chamberlin ff(| to explain this universal behaviour, the diffusive method of Perez-Madrid [l7| . 
CNJ \ and the new method recently proposed by Schwartz and Edwards for non-linear fields. 

In this work, we present an alternative way to obtain the correlation function and a result which is general as long 
as the diffusive processes can be described by a generalized Langevin equation (GLE) . We derive a correlation function 
which is an even function of time. The asymptotic behaviour is similar to those obtained using fractional derivatives. 
However, for short times, our method can lead to richer behaviours. The paper is organized as follows: first, we give 
the outline of the problem of a system governed by a GLE and show the connection to diffusion; we then continue by 
defining in a clear manner the noise, memory and correlation functions. We construct the memory function from the 
fluctuation-dissipation theorem using a coloured noise, from a given density of states. After that, we discuss normal 
and anomalous diffusions. Finally, we analyze the difference between Markovian and non-Markovian behaviours. 

II. GENERALIZED LANGEVIN EQUATIONS AND DIFFUSION 
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We shall start by writing the GLE for an operator A 



dA(t) 



= - f Tl(t-t')A{t')dt' + F(t), (1) 
Jo 



dt 

where F(t) is a stochastic noise subject to the conditions (F(t)) = 0, (F(t)A(0)) = and 

C F (t - t') = (F(t)F(t')) = {A 2 ) eq Ii{t - t'). (2) 

Here, the angular brackets denote an ensemble average. Some correlation functions depend on two times, being 
consequently non-stationary fj, 0, ■ Equation (J2J) is the Kubo fluctuation-dissipation theorem (FDT) [ilj, l2fl | . The 
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FDT can be violated in many slow relaxation processes |2lJ, |22|. Before continuing, it is important to stress the 
connection between the generalized Langevin equation for a given operator A and diffusion. To do this, we define the 
variable 



x{t) = f A{t')dt'. (3) 
Jo 



Now, we can study the asymptotic behaviour of its second moment lim t ^ OQ (a; 2 (i)) ~ t a to characterize the type of 
diffusion presented by the system: for a = 1 , we have normal diffusion, and for a < I (> 1), subdiffusion (superdiffu- 
sion) . 

The memory kernel n(t) indicates that the motion is non-Markovian; when H(t) = 2jS(t), eq. becomes the 
usual Langevin equation. Our main interest is to discuss the behaviour of the correlation function 

_ : A(t)A(0) > 
{) <A(0L4(0)>' 

from which we can describe most of the processes of interest, including relaxation. We use the above conditions for 
the noise to obtain a self-consistent equation for R(t) 



dR{t) 



= - f R{t-t')Ii{t')dt'. (4) 
Jo 



dt 

If we then apply the Laplace transform (denoted by a tilde) to eq. Q , we get 

R(z) = —I—. (5) 
z + U(z) 

From the analysis of this equation, it is possible to obtain plenty of information concerning the asymptotic behaviour 
of the system. In order to make progress, we need to make some assumptions about the origin of the memory. The 
direct way is to connect the random force, F(t), to a thermal bath composed of harmonic oscillators 0, |23, |24| . 
Therefore, for a system in contact with a heat reservoir (canonical), the memory becomes 



U(t) = J p{u>) cos(wt) duo, (6) 

where p is the noise density of states. The memory is clearly even for any noise distribution. We shall use a coloured 
noise given by a generalization of the Debye spectrum 

>-\^) 3 > ifw <^ (7 ) 
, otherwise, 

with uc as a Debye cutoff frequency. The motivation for considering such cases arises from previous studies pil |23| 
in which it was proved that if II(z) cx z v as z — > 0, then the diffusion exponent is given by 

a = l + v. (8) 




III. MEMORY, CORRELATION FUNCTIONS AND NOISE 

In order to obtain R(t), we state that its derivative must vanish at the origin, due to eq. therefore, the 
correlation function can be neither an exponential nor a stretched exponential. The analytical continuation of the 
Laplace transform of an even function is an odd function, and vice- versa. From eq. |J5}, we see that R(z) is odd, 
because n(£) is even (see eq. (©)■ Following the same line of thought, we arrive at the conclusion that R(t) is even. 
Lee [25l | has also shown that both the memory and the correlation function must be even functions for any Hamiltonian 
system. Consequently, we can write 



oo 

R(t) = £ a n t 2n , 

n=0 



(9) 
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FIG. 1: Plot of the correlation function R(t) for normal diffusion (u — 0). Here, 7 = 1 and u>d = 2 and 20, for curves a and b, 
respectively. Curve c is the plot of exp(— 7^;). Note that curves b and c are indistinguishable, except near the origin (see inset). 
Inset: Plot of curves b, c and cos(o;o£), near the origin. 



with ao = -R(O) = 1 and, 

00 

n(<) = ^M 2 ". (10) 

71=0 

We replace those, eqs. and (Q) in eq. (J2J to obtain the following recurrence relation 



2 7 w fl ^ (-l)'[2(n-l-Q]!^ 
°" - -^0! ^ (27+1+^0 a - 1 -'' 

which displays a complex behaviour, where every order depends on all previous ones. This is not a surprise for a 
non-Mar kovian system and it is very useful to have it made explicit in a power series. 

The power series defined by eq. 10 is convergent: it is straightforward to prove by induction that r n — a n /a n+ \ < 0. 
Besides that, it is easy to show that the absolute values of a n decrease for n > hq, for a certain no. Consequently, 
since eq. (0 is an alternating series, and r = linin^oo \r n \ — > 00, it is possible to conclude that it converges |2fJ for 
every value oft. By playing with the parameters v, Wn, and 7, it is possible to generate a large number of functions. 



IV. NORMAL DIFFUSION 



Before we analyze anomalous diffusion, it is very useful to give a more detailed description of normal diffusion. If 
one takes eq. (J7J) with (3 = and ujjy finite, diffusion is normal, despite the noise being coloured. Using this noise in 
eq. ©, we get H(t) = (2j/Tr)s'm(w s t)/t. Applying the Laplace transform yields LT(z) = (2j/ir) arctan(w s /z). It is 
interesting to note that this memory function is realized in electron gas models |27j . We now show that it is possible 
to find at least three time scales for this memory. The first one is the short time scale, given by td = ■ In the 
short time scale, given by t -C t s , the first term in the expansion, eq. @, gives 

R(t) - cos(ujot), (12) 

where uiq = \/n(0)- For this particular case, 11(0) = 2joj s /tt introduces the intermediate time scale r* = ujq 1 oc 
\Jtqtd- In the long time scale, t ^> td, or, equivalently, z <§; cu s , then H(z) ~ 11(0) = 7. Therefore, from eq. JHJ, 
we obtain R(t) ~ exp(— i/ T o)i with To = 7" 1 - This result is the same as the one obtained from the normal Langevin 
equation, which can be obtained from the limiting case u>d — * 00. 

In fig. n we plot the correlation function R(t) as a function of t for normal diffusion. We use 7 = 1, and oju = 2 and 
20 in curves (a), and (6), respectively. Curve (c) is the plot of the function exp(— "ft). It is clear that curve (a) is not 
an exponential; however, in many practical situations it is fine to approximate curve (b) by an exponential for times 
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FIG. 2: Log-log plot of the correlation function R(t) for ujd/^/ = 20 and v = 0.9, 0.5, and 0.1 from top to bottom. Notice the 
transition from a stretched exponential behaviour to a power law. 

larger than tjj. For short times, t <C td, eq. (|12|l is the solution and we disregard the exponential as the solution 
in both cases. The inset shows the correlation function R(t) near the origin for ujd = 20. Together with curve (b) 
and (c), we plot cos(woi). In short, "good exponentials", such as curve (c) cannot exist as a solution in this case. In 
many experimental situations, the ratio Us/l can be very large. Nevertheless, we propose that there will always be a 
short time behaviour which permits distinction from an exponential. Even if the distinction is not accessible to the 
experimental probe, it is very important to bear in mind that the correlation function is even. 

V. ANOMALOUS DIFFUSION 

Now let us return to the general memory. By substituting eq. Q into cq. © with /? ^ and calculating the 
Laplace transform, we get 



n(z) = ^ p (u./z) 

7T 



(13) 



Here. 



V (x) = x-P 



o 1 + 2/ 



(14) 



where a,b) is the Lerch Transcendent 28]. This function has been found most useful in certain quantum statis- 
tics j2^| . This memory function has an asymptotic behaviour of the following form 



hm n(z) ~ <- V, 

2^0 



(15) 



where 



(3, for (3 < 1 
1, otherwise. 



(16) 



This cutoff shows that the exponent a, eq. (JSJl, cannot be larger than 2 (ballistic motion). Recent results [24| confirm 
that cutoff; they obtain v = 1 when (3 = 2. 

By inverse transforming cq. (0), using eq. (|15|l . we get 



R(t) ~ E^i-it/r,) 1 -"), 



(17) 
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which is a useful approximation for t > tjj. Here Eg (x) is the Mittag-Leffler function |13l l30| . This function behaves 
as a stretched exponential for short times and as an inverse power law in the long time regime (see fig. El- Due to the 
asymptotic power-law behaviour 

Ei- u (-(t/T v ) 1 -') ~ [{t/T v ) l -»T{-v)]-\ (18) 

for < v < 1, we obtain an infinite diffusion constant D, which is a sign of superdiffusion |23|. The power law 
exponent 1 — v = 2 — a is in the interval (0, 2) as observed long ago by Friedrich Kohlrausch. For v = 0, the Mittag- 
Leffler function is the usual exponential and tq is the relaxation time. The limiting cases correspond to v = 1 and 
v = — 1; for these values, we have Eq(x) = (1 — x)^ 1 and E^{x) — cosh(y / x), respectively. This last case corresponds 
to the limit of subdiffusion, i.e., a = 0. Since the argument of the function is negative, the Mittag-Leffler function is 
in this case a cosine. For short times, t <C To the first term in the expansion gives a result identical to eq. 1|12|1 . with 
luI = 2ju} D /(l + v)tt. 

In fig. |3 we display the correlation function R(t) as a function of t for anomalous diffusion, obtained from the 
numerical integration of eq. @. We select broadband noise with ujdH — 20. From top to bottom, v — 0.9,0.5 and 
0.1. We notice the transition from stretched exponential to power law behaviour. Note as well that as the behaviour 
becomes more anomalous, i.e. as v increases, the relaxation becomes slower. 



VI. MARKOVIAN VS. NON-MARKOVIAN BEHAVIOUR 



Recently some experimental situations have been found where non-Markovian processes display Markovian be- 
haviour after a long time. In particular, we call attention to the experiments of Merikoski et al. about flame 
evolution [21|. We show now that this transition is possible for normal diffusion. For large t, eq. (@J can be written 



as 



dR(t) 
dt 



-R(t) f Tl(t')dt'. (19) 
Jo 



We can justify this approximation by the following arguments: if t' <C t, then R(t — t') R(t) and for those values 
Il(i') gives the major contribution. As t' increases, say t' f=a t/2, both R and II are small. Finally, for t' w t, 
R(t - t') « 1, but n(i') is very small. 

Notice that this decoupling happens in many situations in physics, otherwise all correlated systems would present 
aging. For normal diffusion, the convergence of the integral of the memory function to a constant value Il(z — > 0) = 7 
makes the decoupling acceptable. In that case, for broadband noise R(t) ~ exp(— jt), as expected (see fig. QJ. 
Consequently, as long as the diffusion is normal, 7 = 11(0) is finite independently of the character of II(t); after a 
time t > 7 _1 all the processes behave as those governed by the usual Langevin equation. For anomalous diffusion, 
the dependence on time of the integral of the memory in eq. I|19fl makes the decoupling fail. In short, the exponential 
behavior in normal diffusion is associated with time invariance while the non-exponential behaviour is associated with 
an aging process which is in essence non-Markovian. 



VII. CONCLUSION 



We discussed the relaxation properties of a dynamical variable A(t) and we proved that neither exponentials nor 
stretched exponentials can give a full description of the relaxation process. However, these functions can be the 
asymptotic behaviour found in many experimental and theoretical works and we found an even function that can 
approach an exponential as we broaden the noise (increase the ratio lod/j). For broad noise, the correlation function 
for anomalous diffusion will be a stretched exponential and finally a power-law. 

Fractional derivatives have been successfully applied to the study of subdiffusive behaviour Those yield 

a = 1 + v, for — 1 < v < 0. The GLE formalism allows v to assume the values in the range pll l23l l32]| — 1 < v < 1, 
which encompass both subdiffusion and superdiffusion. There is a great advantage in our formulation, since the noise 
density of states exists not only in systems governed by a GLE, but in most physical systems. Also, we avoid the levy 
function which has infinite mean square deviation. 

As we have already mentioned, glassy systems seem to be a rich field for studying these phenomena. It would be 
very helpful if the relaxation process could be studied and exponent a determined for those diffusive processes. Other 
related phenomena are chaos synchronization |33l l34| and anomalous reaction rates |35 j. I3& l37l l38j| . which we expect 
to discuss soon. 
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